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QJ • Abstract 
I 

We extend previous analytical studies of the ground-state phase diagram of a one-dimensional Heisenberg spin chain coupled to 
^ . optical phonons, which for increasing spin-lattice coupling undergoes a quantum phase transition from a gap-less to a gaped phase 
with finite lattice dimerisation. We check the analytical results against established four-block and new two-block density matrix 
renormalisation group (DMRG) calculations. Different finite-size scaling behaviour of the spin excitation gaps is found in the 
adiabatic and anti-adiabatic regimes. 
© 2008 Elsevier B.V. All rights reserved. 
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Quantum phase transitions in low-dimensional corre- 
lated systems have attracted considerable attention over 
the last decades. An example is the Peierls instability in 
quasi-one-dimensional spin systems, i.e., the instability of 
a uniform spin chain towards dimerisation induced by the 
interaction with lattice degrees of freedom. Starting in the 
seventies with organic compounds [1], the interest in the 
Peierls instability was renewed with the discovery of a spin- 
Peierls (SP) transition in the inorganic compound CuGeC>3 
in 1993 by Hase et al. [2]. The most significant feature 
distinguishing CuGeOa from other SP-compounds is the 
high frequency uj of the involved optical phonons, which is 
comparable to the magnetic exchange interaction J. 

As an archetypal model for this type of SP system we 
consider the antiferromagnetic Heisenberg chain coupled to 
harmonic Einstein oscillators, 

H = jY^Si-Si+x + uj^lh 

i i 

+ gJ2( b i+ b i)( S i- S i+l- S i- S i-l)> 

i 

where Si denote spin-^ operators at lattice site i, and b\ and 
hi are phonon creation and annihilation operators, respec- 
tively. This model has been studied with a number of ana- 
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lytical and numerical methods, including second-order flow 
equations [3] , the combination of Schrieffer- Wolff transfor- 
mations [4] with a variational ansatz [5], as well as exact 
diagonalisation [5] and a four-block variant of the density 
matrix renormalisation group (DRMG) [6] . All these stud- 
ies agree on the main finding that for finite phonon fre- 
quency the system undergoes a transition from a phase with 
gapless spin excitations to a dimerised phase with massive 
spin excitations only for some finite value of the rescaled 
spin-phonon coupling g = g/u>. 

In the anti-adiabatic limit u> —> oo this critical coupling 
is approaching zero, whereas in the adiabatic limit u> — > 
flow equations and DMRG point towards a finite limiting 
value of g c . This is equivalent to the bare critical coupling g c 
approaching zero linearly with to. To second order in g, the 
Schrieffer- Wolff approach is equivalent to the flow-equation 
result, but when fourth order terms are included the critical 
coupling diverges for decreasing lo. This finding, of course, 
cast doubts on the method and on the decoupling of spins 
and phonons by unitary transformations in general. 

Using improved spin algebra codes in this article we 
therefore reconsider the Schrieffer- Wolff approach for the 
model of Eq. ( 1 ) . In more detail, we try to decouple spin and 
phonon degrees of freedoms by applying a unitary trans- 
formation H = cxp(S)H cxp(— S) that removes interaction 
terms linear in g, 
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Expansion of the effective exchange interactions J n /J up to order 8 
in the electron phonon coupling g. 
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S = gY,(b\-bi){S i -S i+1 -S i -Si- 1 ). (2) 

i 

Averaging the resulting Hamiltonian over the phonon vac- 
uum we obtain an effective spin model with long-ranged 
Heisenberg interactions. Such frustrated spin chains are 
known to be susceptible to dimerisation, and we can use 
the ratio of the next-nearest-neighbour to the nearest- 
neighbour exchange, a e s, as an indicator of the phase 
transition, which for the frustrated spin-^ chain 

H = Y,J(S i -S i+1 +aS i -S i+ 2) (3) 

i 

occurs at a c — 0.241167 [7]. The transformed Hamiltonian 
H can be expanded in a series of iterated commutators, 
H = Efc[S,if] fe /fc!, where [S,H] k+1 = [S,[S,H] k ] and 
[S, H]o = H. For orders k > 2 the evaluation of these com- 
mutators rapidly becomes complicated, and is feasible only 
with efficient computer algebra tools. We were now able to 
push the limit of the expansion from order 4 to order 8. 
Neglecting terms with more than two interacting spins, the 
resulting effective spin Hamiltonian reads 

5 

i n—1 

with J„ = J Y^j=o c i- n ff 2 "* an d the coefficients Cj. n collected 
in Table 1. 

In Figure 1 we show the phase diagram obtained from the 
condition a c s = J2/J1 = a c . The resulting critical coupling 
g c oscillates with increasing order in g. Here the second and 
fourth correspond to the known results of Refs. [3] and [5] , 
respectively, but convergence is achieved only when going 
well beyond that. 

It turns out, however, that the lowest order result shows 
the best agreement with the four-block DMRG data of Bur- 
sill ct al. [6]. On the one hand, this may not seem surprising, 
since the concept of integrating out phonon degrees of free- 
dom usually is appropriate only in the anti-adiabatic limit. 
On the other hand, the smallest energy scale in the prob- 
lem always corresponds to the spin degrees of freedom, and 
therefore an effective spin model should be able to describe 
the phase transition. 

Numerically, the value of g c is obtained from the level 
crossing of the lowest singlet and triplet excitations of the 
full spin-phonon model, Eq. (1), i.e., from the same crite- 
rion applied to get a c of the frustrated spin chain, Eq. (3). 
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Fig. 1. Upper panel: Ground-state phase diagram of the SP model, 
Eq. (1). Lower panel: System size dependence of the level crossing 
at g c in the adiabatic and anti-adiabatic limit. 

In the adiabatic limit this procedure is rather delicate, since 
the appropriate (finite-size) gaps must be smaller than u>, 
requiring large system sizes. In addition, g c has a notice- 
able system size dependence, which is negligible in the anti- 
adiabatic case (see lower panels of Figure 1). To cross-check 
this numerical data we performed new, large-scale parallel 
DMRG calculations [8] using the standard lattice growth 
method of two sites per iteration. For comparable lattice 
sizes we find very good agreement. 

How the unitary transformation scheme can be modified 
to yield improved effective Hamiltonians for the low energy 
physics of spin-Peierls systems will be the subject of future 
studies. 

We thank HLRN and NERSC for providing computa- 
tional resources. 
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